Azimuthal Dependence of Pion Interferometry at the AGS 
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Two-pion correlation functions, measured as a function of 
azimuthal emission angle with respect to the reaction plane, 
provide novel information on the anisotropic shape and ori- 
entation of the pion-emitting zone formed in heavy ion colli- 
sions. We present the first experimental determination of this 
information, for semi-central Au+Au collisions at 2-6 AGeV. 
The source extension perpendicular to the reaction plane is 
greater than the extension in the plane, and tilt of the pion 
source in coordinate space is found to be opposite its tilt in 
momentum space. 

Nuclear collisions at finite impact parameter are in- 
trinsically anisotropic with respect to the azimuthal an- 
gle about the beam direction. For nearly two decades, 
a growing community has studied the details of the 
momentum-space anisotropy of emitted particles (di- 
rected and elliptic flow) at all collision energies [jjj. 
Theoretical studies |^-|7j suggest that coordinate- space 
anisotropics are equally interesting. Recently, a for- 
malism was proposed to use azimuthally-sensitive inten- 
sity interferometry, correlated event-by-event with the 
reaction plane, to extract this information experimen- 
tally j^Q. In this Letter, we present the first experimen- 
tal measurement of the full coordinate-space anisotropics 
of a hot nuclear source, from pions emitted from Au+Au 
collisions at 2-6 AGeV. 

Two-particle interferometry (HBT) is the most di- 
rect probe of the space-time structure of the hot sys- 
tem formed in high energy heavy ion collisions [^],[|, 
and, along with flow, the past decades have also seen 
extensive pion HBT measurements map out the space- 
time geometry and dynamics of heavy ion collisions from 



Bevalac (~1 AGeV) energies @, at the AGS (2-11 
AGeV) 0-0, and to the highest available energy (~160 
AGeV) at the CERN SPS f^ PI W ith the exception of 
some recent preliminary studies 1 12 l^Jl^l, all experimen- 
tal studies to date have assumed cylindrical symmetry of 
the source about the beam axis, strictly valid only for 
6 = 0. 

The E895 Collaboration at the AGS used a large accep- 
tance detector to measure charged particles from Au+Au 
collisions. For every event, the charged particle multiplic- 
ity is used to estimate the the magnitude of the impact 
parameter vector |b|, with an uncertainty ^0.5 fm. For 
all data presented here, the estimated impact parame- 
ter range is |b|=4-8 fm. The direction of b is determined 
from the azimuthally anisotropic momentum distribution 
of protons and light nuclear fragments (not pions), with 
an estimated uncertainty of ~ 20 — 35°. Further details 
are available elsewhere [|l8| , p~9| . 

Experimental details, such as correction for mo- 
mentum resolution and particle identification, of the 
(azimuth-integrated) HBT analysis have been pub- 
lished jll 12 1 . Here, we discuss points particular to the 
azimuthally-sensitive analysis. 

As in "standard" HBT analyses, correlation functions 
are constructed by dividing the two-pion yield as a func- 
tion of relative momentum q = pi — P2, by a background 
generated by mixing pions of the current event with those 
of previous events [ pp[ . As is well-known, single-particle 
acceptance/efficiency and phase-space effects cancel out 
with this "event-mixing" technique. However, there is a 
special consideration here. In the present analysis, we 
generate correlation functions with selections on the pair 



1 



angle with respect to the reaction plane <f> — ^(Kj_,b), 
where Kj_ = (pi + P2)± is the total momentum of the 
pair perpendicular to the beam H). Thus, we only mix 
events which have similar (within 5°) directions of recon- 
structed b. 

The correlation function was binned according to q- 
components in the Bertsch-Pratt ( "out-side- long" ) de- 
composition as measured in the Au+Au cm. frame [ p2[ . 
In this Letter, we use the subscripts "o, s, 1" to stand 
for "out, side, long." Here, qi is the component parallel 
to the beam, q a is the component parallel to K, and q s 
is perpendicular to qi and q Q . Particles in a pair were 
ordered such that qi > 0; the signs of the other two com- 
ponents were retained. 

Figure |l| shows two-dimensional projections of the cor- 
relation function for 4 AGeV Au+Au semicentral colli- 
sions, measured in eight 45°-wide bins in <f>. The average 
rapidity and pr of pions in low-g pairs was 1.2 (midra- 
pidity) and 110 MeV/c, respectively, and did not vary 
with 0-bin. 

Especially in the q ~qi and q s —qi projections, the cor- 
relation functions display a distinct "tilt," which evolves 
with the emission angle <j). As discussed below, this 
tilt arises from fundamental geometry 0, in contrast 
to a similar tilt observed JLsf in some q Q — qi projec- 
tions measured far from midrapidity for azimuthally- 
integrated HBT, which arises from dynamic (longitudi- 
nal flow) effects p3fl . Aside from these dynamical effects, 
there are no 2-dimcnsional tilts in azimuthally-intcgrated 
HBT, and only one-dimensional projections are typically 
shown JE^ |. 

Using a maximum-likelihood technique, we fit the cor- 
relation functions for each of the <fc bins with the standard 
Gaussian parameterization pi 



C(q,. 



1 + A(</>) exp 



i } j=o,s,l 



(1) 



In contrast to azimuthally-integrated analyses, all six ra- 
dius parameters are relevant here . 

Two-dimensional projections of the fits, weighted ac- 
cording to the mixed-event background, are shown as 
contours in Figure ^. Since the relative sign of any 
component of q = pi — P2 is arbitrary (while the rela- 
tive sign between components is not arbitrary), the two- 
dimensional projections in Figure 1 have been reflected 
for clarity. However, the fits were performed in the full 
three dimensions, with each pion pair being weighted 
only once. 

Although the full structure of the correlation function 
is only visible in multi-dimensional projections such as 
those in Figure ^, one-dimensional projections are com- 
monly shown for azimuthally-integrated HBT analyses, 
to provide a feeling for the quality of the data and fit. 
For reference, one-dimensional projections of the corre- 
lation function and fits, for one <j> slice, are shown in 
Figure ||. 



The seven fit parameters Rij and A are plotted as a 
function of (j) in Figure |. As is well-known [§,§, the 
HBT "radii" parameters are in general not strictly the 
radii of the emitting source, but have a geometrical in- 
terpretation as lengths of homogeneity of the emitting 
region. For realistic sources, various interpretations of 
A have been proposed (see e.g. p^ , |25[ |). At AGS ener- 
gies, the observed strong energy dependence of A was 
well reproduced O] by detailed simulations using the 
RQMD (v2.3) transport model |^(|; there, the decrease 
of A with energy was attributed to increased production 
of long-lived 7r _ -emitting particles at higher energy (see 
also f27|] ). In this scenario then, the HBT radii charac- 
terize a Gaussian "core" of directly-produced pions, with 
A < 1 resulting from a well-separated "halo" of pions 
emitted from long-lived particles |28| . In this framework, 
we now discuss the ^-dependence of the HBT parame- 
ters, and how to extract the underlying geometry of the 
"core" of the source. 

The </>-indepdence of A indicates that the fraction of 
7r~ from long-lived resonances is independent of emission 
angle with respect to the reaction plane. This seems 
reasonable, given the relatively weak anisotropic flow at 
these energies Jl8|] , as well as the randomizing effect on 
the pion from decay of the parent. 

R and R s show significant equal and opposite second- 
order oscillations in </>, consistent with the simple picture 
of a transverse profile of the pion-cmitting region that re- 
flects the "almond-shaped" overlap region (with a larger 
spatial extent perpendicular to b than parallel to b) be- 
tween the target and projectile spheres ]3|,|l2| • 

Note that the "cross-term" radii (-R 2 ;, R^, R 2 sl ) quan- 
tify the tilt observed in the 2-dimensional correlation 
function projections of each 0-bin of Figure |l|. R 2 ol 
and R 2 sl display equal-magnitude first-order oscillations, 
consistent with pion emission from an ellipsoidal source 
which is tilted in coordinate space in the reaction plane, 
and away from the beam axis ||. 

In general, the six HBT radii are related to the spa- 
tiotemporal structure of the emitting source via the equa- 
tions nu 



R„ 



= Sii sin 2 (f> + S22 cos 2 <; 
= S11 cos 2 (j) + S22 sin 2 <; 

-2/3 ± S Q1 cos(f>-2P_ 
= S33 — 2(3iSq3 
, = 5i 2 cos2</>+ \ (S22 
+/3j_<S'oi sin0 - (3 



— S12 sin 2</> , 
+ S12 sin 2(f> 
„S 02 sin0 + P±S 



f3?Soo ■ 



- Sn) sin 2<j) 
S 2 cos (j) , 

Rli = (Sw-ptSoJcaat-PxSaa 

+ (5*23 ~ /3i S a2 ) sin (f> + f3iP±S Q0 , 
Rli = (S23 ~ AS02) cos cj> - (S 13 - AS01) sin <t> . 



(2) 



where fix and j3i are the average pair velocities in the 
transverse and longitudinal direction, and the spatial cor- 
relation tensor is given by 
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x„, = 0,1, 2, 3) (3) 



where the brackets ( ) indicate an average over the emit- 
ting source. In Equations @, fi = 3 refers to the beam 
direction, and /i — 1 is parallel to b. 

Simultaneously fitting the six i?? ■ {(f) of Figure |^ with 
Equations^ (treating the as 10 (^-independent fit pa- 
rameters) allows extraction of the full spatial correlation 
tensor 



/ 13.4 ±3.5 
0.7 ±0.9 
0.2 ±1.0 
\ 0.0 ±1.1 



0.7 ±0.9 
21.3 ±0.9 
-0.1 ±0.6 

4.5 ±0.6 



0.2 ± 1.0 
-0.1 ±0.6 
25.5 ± 1.0 
-0.3 ±0.6 



0.0 ± 1.1 ^ 
4.5 ±0.6 
-0.3 ±0.6 
22.8 ±0.9 J 



where all units are in fm 2 and c = 1 . The fits are shown 
by solid lines in Figure |3|. 

The only significantly non-vanishing tensor elements 
are the diagonal ones and 6*13. This fact, along with 
observations of identical oscillation amplitudes in i? 2 and 
R 2 , and in R 2 t and R 2 sl , support the assumption of (f>- 
independence of S^ IV in the fit |||(| . 

S13 and (622 — S11) quantify, respectively, first and sec- 
ond harmonic oscillations of the HBT radii with respect 
to the measured reaction plane. However, due to finite 
particle multiplicity, the measured reaction plane differs 
statistically from the true one by some angle A</> ]l8|,^9| . 
This reaction plane dispersion results in an apparent re- 
duction in both 5*13 and (S22 — Sn) (but does not af- 
fect quantities Sqq, S33, and (Sn + 822))- According 

to MMM 



Si3, m = S13 ■ (cos(A0)) (4) 
(S22 - S n ) • <cos(2A0)>, (5) 



(5*22 — Sn 



where 5i3 iTO is the observed value, and S13 is the "true" 
value. The correction factors obtained by following the 
procedure outlined in Refs. @JT|,||,|| are listed in Ta- 
ble |. Hence, for the 4 AGeV data, the reaction plane- 
dispersion corrected values are Sn = 19.8 ± 1.2, fm 2 , 
S 22 = 27.0 ± 1.4 fm 2 , and S13 = 5.2 ± 0.7 fm 2 . 

The spatial tilt of the emission ellipsoid in coordinate 
space is given by || 



1 

— tan 
2 



2S 



13 



S: 



Sn 



37° ± 4°. 



(6) 



Rotating the spatial correlation tensor (corrected for 
reaction-plane dispersions) by this angle returns the di- 
agonal tensor 



( 13.4 ±3.5 
0.5 ± 1.0 
0.2 ± 1.0 
\ 0.4 ±1.0 



h)-S- Ry{ 

0.5 ± 1.0 
15.9 ± 1.0 
0.1 ±0.6 
0.0 ±0.7 



0.2 ± 1.0 
0.1 ±0.6 
27.0 ± 1.4 
-0.3 ±0.6 



0.4 ± 1.0 \ 
0.0 ±0.9 
-0.3 ±0.6 
26.8 ±0.8 J 



whose diagonal elements are the squared lengths of ho- 
mogeneity {Spy = a 2 ) in the t, x, y, z directions. 

Similar results are obtained at other energies. Fig- 
ures [| and |5| s how the HBT parameters obtained at 2 
and 6 AGeV |3l],|3^]. Especially at the highest energy, de- 
graded reaction-plane resolution significantly reduces the 
second-order oscillations in the transverse radii. Consid- 
ered in isolation, R 2 in Figure |B| may be fit by a constant 
(^-independent) value as well as anything else. How- 
ever, the simultaneous fit of all R 2 (4>) allows relatively 
clean extraction of the source shape parameters. Reac- 
tion plane dispersion corrections are applied to all data, 
but the effects are small. For the worst case (6 AGeV 
dataset), the corrections result in a decrease (increase) 
in a x (a y ) of 0.3 fm, and 5° reduction in 8 S ; a z and at 
are unchanged. 

Table H summarizes the inferred pion source shapes for 
our three energies. At all energies, our results indicate 
a pion freezeout distribution as an ellipsoid whose major 
axis in the reaction plane is tilted with respect to the 
beam in the positive direction (i.e. in the direction of b), 
and whose transverse axis perpendicular to the reaction 
plane is longer than the axis in the reaction plane (the 
"almond" shape referred to above) . The extension in the 
temporal direction is consistent with observations at high 
energy 

The observed transverse "almond" shape is reminis- 
cent of the entrance channel geometry. In the simplest 
picture, pions are emitted from the overlap region of 
the two spherical Au nuclei. For impact parameter b = 
4 fm (8 fm), the spatial RMS of the overlap region is 
2.2 fm (1.3 fm) parallel to b, and 2.9 fm (2.4 fm) per- 
pendicular [[33]]. While the linear scale of the freezeout 
distribution is roughly twice these estimates (see Ta- 
ble n, indicating significant expansion, the aspect ratio 
cr y /a x m 1.3 — 1.4, is in the range (1.3 — 1.9) expected 
from this naive picture. 

On the other hand, the large positive tilt is a geomet- 
ric feature of the collision dynamics ||. Studies with 
the RQMD model indicate that most of the 1ow-pt pions 
arise from A decay, so it is not surprising that the tilted 
freeze-out distribution resembles the baryonic distribu- 
tion in coordinate space calculated in transport codes [|J . 

Experimental access to this level of geometric detail 
on the freezeout distribution is unprecedented, and rep- 
resents an exciting new opportunity to study the dynam- 
ical response of hot nuclear matter to compression. A de- 
tailed theoretical discussion is beyond the scope of this 
Letter, but we note that an identical analysis on pions 
generated by the RQMD model displays considerable sen- 
sitivity on the dynamical effect of the nuclear meanficld; 
RQMD values are listed in Table [El]. Although the spa- 
tial scale is underpredicted at the lower energies and the 
temporal scale overpredicted (noted already for central 
collisions |ll|), qualitatively, the model reproduces the 
"almond" shape and large positive tilt angles remark- 
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ably well. Since the model better describes proton di- 
rected flow when the mean field is included in the cal- 
culation Q , it is interesting to note that the tilt angles 
(the spatial counterpart of proton directed flow 0) re- 
produce observation better when the mean field is off. 

We stress that these coordinate-space anisotropies rep- 
resent new information, independent of momentum-space 
anisotropies (directed and elliptical flow). Momentum- 
space tilt angles (flow angles Q) at these energies are 
only a few degrees 19 , and, indeed, the directed flow 
(momentum-space tilt) of pions at these energies is in 
the negative direction p5| , p6t p[, opposite the coordinate- 
space tilt 8 S . Experimental information on the interplay 
between coordinate and momentum space anisotropies 
should help resolve theoretical issues, such as the coex- 
istence of flow and antiflow components at these ener- 
gies g. 

Further studies of the elliptical transverse shape of the 
source at higher energy should prove quite interesting. It 
is known, for example, that the elliptic flow of nuclear 
matter (protons and pions) in momentum space changes 
sign at AGS energies [ ^6|Jl^ , |37| , |3^ | , from negative (more 
momentum out-of-plane) at low collision energy to pos- 
itive. The coordinate-space analogue of negative ellip- 
tic flow is the "almond" shape we observe. The degree 
to which this shape follows the momentum space and 
evolves to an in-plane shape at high energies, may provide 
valuable information on the detailed nature and cause of 
elliptic flow RQMD predicts Q that the "almond" 
shape is retained even at RHIC energies. 

In conclusion, we have presented the first full measure- 
ment of the azimuthal dependence of pion interferome- 
try. For semi-peripheral Au+Au collisions at 2-6 AGeV, 
the spatial correlation tensor, S^, extracted from the 4>- 
dependences of the HBT radius parameters, reveals an 
ellipsoidal pion emission region with an "almond" trans- 
verse profile and which is strongly tilted in the reac- 
tion plane away from the beam axis. Consistency re- 
lations indicate that the extracted geometry is negligi- 
bly affected by possible transverse space-momentum cor- 
relations. The RQMD transport model reproduces the 
qualitative features of the data quite well, and reveals a 
dependence of the coordinate-space anisotropies on the 
action of the nuclear meanfield. Through application of 
this new type of analysis at other energies and careful 
theoretical comparisons, it is hoped that fresh insight on 
the dynamics of non-central heavy ion collisions may be 
gained. 
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FIG. 1. Two-dimensional projections in q — q s (left col- 
umn), q — qi (center column), and q a — qi (right column), 
of the correlation functions measured for midrapidity pions 
from semi-central Au+Au collisions at 4 AGeV. For each 
projection, the unplotted component is integrated over ±30 
MeV/c. Indicated by the labels on the right, projections are 
shown for emission angles with respect to the reaction plane 
(f> = 45° ± 22.5° (top row) to tj> = 360° ± 22.5° (bottom row). 
Solid lines show projections of fits with Equation [lj 
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FIG. 2. One-dimensional projections, in q , q s , and qi, of 
the correlation function for (f> — 45° ± 22.5° (top row of Fig- 
ure |l|). For each projection, the unplotted q components are 
integrated over ±30 MeV/c. Note that in the present analy- 
sis, qi is defined positive, while q and q s have a meaningful 
signs. Solid lines show projections of fits with Equation hi 
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FIG. 3. Filled stars show the fit parameters from Equa- 
tion 111 resulting from fits to the correlation functions of Fig- 
ure \^ The values at <j> — 0° are redisplayed as open stars 
at 4> = 360°. The line in the A panel represents the average 
value of A. Lines in the other panels represent the fit to the 
HBT radii (stars) with the Equation The vertical scale is 
linear in all panels. 
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FIG. 4. Same as Figure Bl but for 2 AGeV collisions. 
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FIG. 5. Same as Figure til but for 6 AGeV collisions. 



E (AGeV) 


(cos(A0)) 


(cos(2A<j!>)) 


2 


0.940 


0.787 


4 


0.853 


0.584 


6 


0.724 


0.384 



TABLE I. First and second harmonic reaction plane dis- 
persion correction factors for the three datasets discussed. 
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4 Data 


3.7 ± 0.5 


4.0 ± 0.1 


5.2 ± 0.1 


5.2 ± 0.1 


37 ± 4 


[2] 


RQMD cs 


6.4 ± 0.2 


3.3 ± 0.1 


4.2 ± 0.1 


4.7 ± 0.1 


33 ± 3 


RQMD mf 


4.6 ± 0.3 


3.5 ± 0.1 


4.7 ± 0.1 


4.7 ± 0.1 


45 ± 3 


[3] 


6 Data 


3.8 ± 0.5 


3.5 ± 0.2 


4.8 ± 0.2 


4.7 ± 0.1 


33 ± 6 


RQMD cs 


5.9 ± 0.2 


3.5 ± 0.1 


4.3 ± 0.1 


4.9 ± 0.1 


28 ± 3 


[4] 


RQMD mf 


5.0 ± 0.3 


3.6 ± 0.1 


4.7 ± 0.1 


4.5 ± 0.1 


48 ± 5 



TABLE II. Lengths of homogeneity cr M and tilt angle at 
each collision energy. In addition to experimental results, 
predictions of the RQMD model with meanfield off (cs) and 
on (mf) are given. 
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